Possible alterations of the gravitational field in a superconductor 
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In this paper I calculate the possible alteration of the gravitational field in a superconductor 
by using the time-dependent Ginzburg-Landau equations (TDGL). I compare the behaviour of a 
high-T c superconductor (HTCS) like YBa2Cu307 (YBCO) with a classical low-T c superconductor 
(LTCS) like Pb. Finally, I discuss what values of the parameters characterizing a superconductor 
can enhance the reduction of gravitational field. 
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There is no doubt that the interplay between gravitational field and superconductivity is a very intriguing field of 
research, whose theoretical study has been involving many researchers for a long time[[j], 0, ^, ^, ^). Eight years ago, 
E. Podkletnov and R. Nieminem declared the achievement of an experimental evidence for a gravitational shielding 
due to a rotating high-T c superconductor. After their announcement, other groups tried to repeat the experiment 
but they obtain controversial results j|, @j H) > so ^ na * a * t ne present moment the question is still open. 

In 1996, G. Modanese interpreted the results by Podkletnov and Nieminem in the frame of quantum theory of 
General Relativity || but the complexity of the formalism he used makes it very difficult to extract quantitative 
predictions. 

In a very recent paper, M. Agop JlC| ] and collaborators wrote generalized Maxwell equations that simultaneously 
treat weak gravitational and electromagnetic fields. In the weak field approximation, the Einstein equations fnl , 

- -R„^ = Z^Lt^ 
2 y c 4 

lead to the following equations JT2[ , formally similar to Maxwell's: 

V • E 9 = -AnGp g (1) 

V • B g = (2) 

VxE 9 = -^ (3) 
„ „ 4ttG. 1 9E„ 

VxB -- J » + ?¥ (4) 

where E 9 and B 9 are the gravitoelectric and gravitomagnetic field respectively, j 9 is the mass current density vector 
such that j g = vp g , v is the velocity, p g is the mass density. Obviously G is the Newton's constant and c is the speed 
of light in vacuum. As in the electromagnetic case, it is possible to define a gravitational permitivity e g = 1 /AitG 
and a gravitational permeability p g = AttG/c 2 of the vacuum. For example on the surface of the earth E 9 is simply 



the Newtonian gravitational acceleration and B g is related to angular momentum interactions ||10|, 12 . Then, they 
defined generalized electric field, magnetic field, scalar and vector potentials containing both an electromagnetic and 
a gravitational term, in the following way: E = E e +^E S ; B = B e +^B g ; <j) = (j> e +^(j> g and A = A e +^A g where 
m and e are the electronic mass and charge and the subscripts e and g mean 'electromagnetic' and 'gravitational' 
respectively. The generalized Maxwell equations then become [[lO] : 

V-E= (-4nG+ p (5) 



V- B = 



(G) 



2 



V x E = — — 



dB 

dt 



(7) 



VXB=I — + "» 



(8) 



where the relations: p g = ^p and j g = have been used and eo and /io are the electric permitivity and magnetic 
permeability in the vacuum. They also wrote the two generalized London equations |10(| 



E 
B 



(1/p) (dj/dt) 

(-Vp)vxj 



(9) 



and so they could define the generalized penetration depth 

A = 



AgA e 



A 2 -A 2 



(10) 



where A e = [m/(yLt e 2 n)] 



1/2 



A„ 



1/2 

[c 2 / (ATrGmn)] , n is the density of superelectrons and A ff /A e 



10 



21 



For simplicity, I will study the case of an isotropic superconductor, in the gravitational field of the earth, in the 
absence of an electromagnetic field, thus taking E e = and B e = 0. B s in the solar system is very small [|l6| therefore, 
E =^E g and B = 0. Moreover the gravitational effects of j g are insignificant (j g is related to gravitational effect of 



the superconductor) and so I assume j s = 0. I have also ■ 



--^<t> g and A=f A g . 



This situation isn't analogous of the 



Meissner effect but, rather, to the case of a superconductor in a electric field. Since the gravitoelectric fied is formally 
analogous to electric field I will use the time-dependent Ginzburg-Landau equations (TDGL) Jl3|, 14, which, in 
the Coulomb gauge V • A = 0, are written in the form: 



2mD 



i + ^rA i>-aii + b i, + ^-{ihV + —A) 2 t/j = 
at n I 2m c 



(11) 



-47R7 . 1 OA . 4-7T 

VxVxA-VxH = (-^r- + V<£) H 

c c ot c 



eh 



4e 2 9 ' 
— H 2 A 



(12) 



where D is the diffusion coefficient, a is the conductivity in the normal phase, H is the applied field, a(T) = a (T — T c ) 
and b(T) = b(T c ) where ao and b are the positive constants and T c is the critical temperature of the superconductor. 
The boundary and initial conditions are: 



(ihVip + (2e/c)AV0 • n = 

V x A • n = H • n } on dQ. X (0, t) 
A- n =0 



ip(x,0) = ip (x) 
A(x,0)=A (x) 



n 



(13) 



where dfl is the boundary of a smooth and simply connected domain fl in R n . In order to write equations [Ll|,|12| in a 
dimensionless form, the following quantities can be defined: 



* 2 (T) = H C (T) = J^ml 



h/2e 



2^\{T)t(T) 



(14) 



^2m\a(T)\ 



A(T) = 



bmc 2 



47r|a(T)|e 2 



K = A(T)/£(T), r(T) = X 2 (T)/D, V = A7raD/(e c 2 ) 

where A(T), £(T) and H C (T) are the penetration depth, the coherence length and the thermodynamic field. 
The dimensionless quantities are then: 



(15) 
(16) 



x' = x/X, t' = t/r, V' = ip/V 



(17) 
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A' = Ak/(V2H c \), (j)' = (j)K/{V2H c D), H' = Hn/(V2H C 



(18) 



Inserting the eq. [l7 18 in cqs. |TT|,^2| and dropping the prime gives the dimensionless TDGL equations [[l3) in a 
bounded, smooth and simply connected domain SI in R": 



(hp 
dt 



^- + ifiip + k 2 {\iP\ 2 - 1)V> + (iV + A) 2 V = 



rj{^ + Vcj>) + ^i(ip*Vip - ipVip*) + \ipf A + VxVxA-VxH = 



The boundary and initial conditions ( |l3| ) become, in the dimensionless form 

on dtt x (0, t) 



(iS/ip + Aip) • n = 
V x A • n = H ■ n 
A-n=0 



ip(x,0) = i(> (x) 
A(x,0)=A (x) 



on £1 



(19) 
(20) 

(21) 



Our superconductor is immersed in the gravitational field of the earth which is very weak and approximately 
constant. So <p = —g*x where g* = \(T) unrig / (-\/2eH c (T)D) <C 1 and g is the gravity acceleration. The corrections 
to <p in the superconductor are of the second order in g* and therefore they aren't considered here. 

Now I search for a solution of the form: 



ip(x,t) = ipo(x,t) +g*f(x,t) 
A(x,t) = + g*P{x,t) 
<p(x) = —g*x 



(22) 



At order zero in g*, eq.(|l9|) gives: 

dip (x,t) 



dt 



with the conditions: 



V(|VoOM)l -WOM) 

V>oOM = 0) = 
ip {x = 0,i) = 
1>o(x = L,t) = 



d 2 ip {x,t) 
dx 2 



(23) 



(24) 



where L is the length of the superconductor and t — is the instant when the material undergoes the transition to 
the superconducting state. 

The static classical solution of eq. ^3| is: 

(25) 



(26) 



->po(x,t) = ipa(x) — jtanh 
At the first order in g* one obtains from equation ( |l9| ) 



kx/V2 


— tanh 


k(x- L)/V2 


— tanh 


kL/V2 


| / tanh 


kL/V2 



d>y{x,t) d 2 j(x,t) 



dt 



dx 2 



+ k 2 (3 \ip (x) | - l)j(x,t) = ixip (x) 



with the conditions: 



7(x,i = 0) = 
j{x = Q,t) = 
7(jc =L,t) = 



(27) 



The equation at order one for the vector potential is 



+ lMx)l 2 j9 ( a . j t) + j( Xj t )- v = o 



(28) 



with the constraint 



P(x,t = 0) = 



(29) 
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Note that the second-order spatial derivative of (3 does not appear in eq.(ps|). This is due to the fact that, in one 
dimension, V 2 A — • A and therefore, in the Coulomb gauge, V x V x A = V(V ■ A) — V 2 A = 0. The quantity 
J(x, t) which appears in eq. is given by: 



The solution of eq. Eg] is 



J(x,t) 



d d 



(30) 



0(x,t) 



\M*)\ d 



1 - exp(- \ip (x)\ 2 t/Ti) 



exp(- \ipp(x)\ t/rj) 
V 



dtJ(x,t)exp(\ip (x)\ t/rj) (31) 



Now, we have an expression for ipQ(x,t) (eq. Eq), and an expression for 0(x,t) as a function of j(x,t) (eq. pl| ). 
Therefore, we can insert these expressions in eq. |22j and obtain both i/)(x, t) and A(x, t) as functions of j(x, t). Using 
the relation: 



E„ 



we finally find the gravitoelectric field E g in the superconductor: 



E fl (s,t) 
9* 



1 - exp(- \ijJo{x)\ t/ri)-—l 



/ dtJ(x,t)e-xp{\ip (x)\ 2 t/rj) 
Jo 



(32) 



(33) 



From this formula we can see that for maximizing the effect of the reduction of the gravitational field in a supercon- 
ductor it is necessary to reduce rj and to have large spatial derivatives of tpo(x) and j(x, t). The condition for having 
a small value of rj is that the superconductor has a large normal-state resistivity and a small diffusion coefficient 
D ~ vfI/3 (where vp is the Fermi velocity, which is small in HTCS, and I is the mean free path). Therefore, the 
effect is enhanced in 'bad' samples with impurities, not in single crystals. 

From the experimental viewpoint, the greater are the length and time scales over which there is a variation of E g , 
the easier is the observation of this effect. Actually, we started from non-dimensional equations and therefore the 
length and time scales are determined by A(T) and t = X 2 (T)/D, which should therefore be as large as possible. In 



this sense, some new materials with very large A(T) 
as clearly seen in eq. |33|, the relaxation time is inverse 
as possible, and this implies that k is small (see eq 



18j could be interesting for the study of this effect. Moreover, 
y proportional to |-0o(£)| 2 - ^ s a result, ipo{x) must be as small 
Then, A(T) and £(T) must be both large. 
Up to now we have dealt with the expression of ]3(x,t) as a function of j(x,t). Actually, to obtain an explicit 
expression for E g we have to solve the equation for 7(2;, i) (eq. ^6|) . This is a difficult task which can be undertaken 
only in a numerical way. Nevertheless, if one puts tpo(x,t) ~ 1 (good approximation in the case of YBa2Cu30y 
(YBCO), in which re=94.4), one can find the simple approximate solution: 



+00 



f(x, t) — £70(2;) + i A m sin(nirx / L) exp [— (m 2 ir 2 /L 2 + 2n 2 )i\ 



(34) 



where 



7o(x) = x/(2k 2 ) ■ jl - sinh(V2K(L - x)) + swh(V2~K,x) / sinh(V2fi;L)| 



(35) 



If L 

A m = — I 70(1) sin(TO7ra;/L)da; = 

L J ZmirK z 



(-iy 



[A lm + A 2m ] 



(V2KL/(rmr)) +1 



(36) 



and 



Ai m = ran 



1 - (-l) m cosh(\/2«;L) 



(-l) m V2~KLsmh(V2KL) 



(v/2«L/(m7r)) +1 



(37) 
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1 - cosh(\/2KL)l J 2V2k£ [(-l) m cosh(\/2Ki) - l] 
smh(V2nL) \ ) ( TO7r ) (^^/(mTr)) 2 + 



1 



m7r(-l) m sinh(V2KL 



(38) 



By inserting eq. 35 in eqs. [30|j3l| and taking into account eq. ^5], the gravitoelectric field E 5 becomes: 

(x)C m (x,t) 



E g (x,t) = g* 



l-exp(- ^(xft/rfiil-^ 



+ OC 

£ 

m— 1 



where 



J (x) 



1 

2^2 



d d 
4>o{x)—j (x) - j {x)—ip {x) 



(39) 



(40) 



B m {x) = (mir / L)tpo(x) cos(rmrx/L) 



s'm(mirx/L) — i/jq (x) , 



(41) 



(m 2 Tr 2 /L 2 + 2k 2 ) exp [-(m 2 ir 2 /L 2 4 


-2k 2 )^] - (|^o(^)| 2 /^) exp(-|Vo(x)| 2 t/77) 


(\M*)\ 2 /v) 


- (m 2 7T 2 /L 2 + 2k 2 ) 



C m {x, t) 

Note that making the very drastic approximation 

7(x) ~ ix/(2n 2 ) 

leads to the apparently draft result 



E g (x,t)=g* 



l-exp(-\Mx)\ t/v) 1 



^oo(^) 



where 



Joo(x) = i 



■00 (^) - a; 



9Vo(a;) 



dx 



(42) 



(43) 



(44) 



(45) 



In spite of its crudeness, this approximate solution (eq. in the case of YBCO gives the same results of the exact 
solution (eq. |3^). Moreover, nothing changes significantly if one neglects the finite size of the superconductor and uses 
ipo{x) = tanh (nx/\/2) instead of eq. In the case of YBCO the variation of the gravitoelectric field E 9 in time and 
space is shown in Figure la and lb. It is easily seen that this effect is almost independent on the spatial coordinate. 

The results in the case of Pb are reported in Figure 2a and 2b, which clearly show that, due to the very small 
value of k, the reduction is greater near the surface. In this case, moreover, some approximations made in the case 
of YBCO are no longer allowed. For example, for small values of L the simplified relation ( [43| ) is not valid. When k 
is small, in fact, the length L plays an important role. In particular, if L is small the effect is remarkably enhanced, 
as shown in Figure 3. In the same condition, a maximum of the effect (and therefore a minimum of E ff ) can occur at 
t 7^ 0, as can be seen in the same Figure. 

In conclusion, for YBCO the shielding effect decays with a relaxation time i sur f ~ 0.25 r = 8.5 10 -11 s near the 
surface (x = 0.01A = 3.3 • 10~ 9 m) and tj nt — O.lr = 3.4 • 10~ n s in the interior of the sample (x = A = 3.3 • 10 _7 m). 



In the case of Pb, the same quantities take the values i sur f ~ 3 
tint. ^ 5 • 10 5 t = 3.1 • 10~ 9 s (x = 100A = 7.8 • lO^m) with L > A. 



10 't = 1. 



io- 



(x = A = 7- 



Table 1 reports the values of the parameters of YBCO and Pb, calculated at a temperature T such that (T 
is the same in the two materials. 



10~ 8 m) and 
T c )/T c 
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Table 1 





YBCO 


Pb 


T c 


89 K 


7.2 K 


T 


77 K 


6.3 K 


A(T) 


3.3-1CT 7 m 


7.8-10" 8 m 


Z(T) 


3.6-10~ 9 m 


1.7-10~ 7 m 


a- 1 


4-10~ 7 Qm(T = 90K) 


2.5-10~ 9 Qm(T = 15K) 


H C (T) 


0.2 T 


0.018 T 


r(T) 


3.4-10- 10 s 


6.1-10- 15 s 


D 


3.2-10" 4 m 2 /s 


1 m 2 /s 


1 


6-10~ 9 m 


1.7-10 -6 m 


v F 


1.6-10 5 m/s 


1.83-10 6 m/s 


K 


94.4 


4.8-10" 1 


V 


1.27-10" 2 


6.6-10 3 



Tables 2. a and 2.b summarize the variation of the fundamental quantities with the temperature [|lj 

Table 2.a Table 2.b 



YBCO 


A 


T 


9* 


Pb 


A 


T 


9* 


T = OK 


1.7-10~ 7 m 


9.03-10- n s 


2.6-10~ 12 


T = OK 


3.90-10- 8 m 


1.5-10" 15 s 


MO" 17 


T = 70K 


2.6-10~ 7 m 


2.1-10- 10 s 


9.8-10" 12 


T = 4.20K 


4.3-10~ 8 m 


1.8-10" 15 s 


1.4-10" 17 


T = 77K 


3.3-10~ 7 m 


3.4-10- 10 s 


2-10" 11 


T = 6.26K 


7.8-10" 8 m 


6.1-10" 15 s 


8.2-10" 17 


T = 87K 


8-10~ 7 m 


2-10- 9 s 


2.8-10" 7 


T = 7.10K 


2.3-10" 7 m 


5.3-10~ 14 s 


2.2-10- 15 



It is clearly seen that A and r grow with the temperature, so that one could think that the effect is maximum 
when the temperature is very close to T c . However, this is true only for low-T c superconductors because in high-T c 
superconductors (HTSC) fluctuations are of primary importance for some Kelvin around T c . The presence of these 
opposite contributions makes it possible that a temperature r max < T c exist, at which the effect is maximum. 

In all cases, the time constant Ti nt is very small, and this makes the experimental observation rather difficult. Here 
I suggest to use pulsed magnetic fields to destroy and restore the superconductivity within a time interval of the order 
of T int . 

The main conclusion of this work is that the reduction of the gravitational field in a superconductor, if it exists, is 
a transient phenomenon and depends strongly on the parameters that characterize the superconductor. 

Note that in this paper I have used a very simplified model. For a more realistic description, one should take into 
account some features of real superconductors, for example: 

1. The symmetry of the order parameter, which in HTCS can be different from a pure s-wave [ p^[ ; 

2. The fact that the relaxation constant r\ can be complex p(| ; 

3. The high anisotropy and layered structure of HTCS 

4. The effect of superconducting fluctuations, which is very large in HTCS |p2| . 

time and space derivative of the order parameter are present, I suggest 



Finally, since in the general solution (eq. 
that this effect could be enhanced: 

1. by the presence of impurities 

2. by using quickly variable (pulsed) magnetic fields ]24j : 

3. by making the superconductor quickly rotate [^5| ; 

4. by using constant or time-dependent electric fields [ p6| ; 

5. by manufacturing a superconductor made of layers of different materials, or different phases (with different T c ) 
of the same material 27 . 
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FIGURE CAPTIONS 



Fig. 1 (a) The gravitational field E g /g* as a function of the normalized time and space for YBCO at T — 77 K; (b) 
the gravitational field as a function of the normalized time for increasing values of x: x — 5 • 10~ 3 A, x = 10~ 2 A and 
x = X. 



Fig. 2 (a) The gravitational field E s /g* as a function of the normalized time and space for Pb at T — 6.3 K; (b) 
the gravitational field as a function of the normalized time for increasing values of x: x = X, x = 2A and x = 10A. 



Fig. 3 The gravitational field E g /g* as a function of the normalized time in the case of Pb, with L=Q, 8 and 10 and 
x=4. The maximum of the effect is evident. 




FIG. 1: G. A. Ummarino, Possible alterations of the gravitational field in a superconductor 




FIG. 2: G. A. Ummarino, Possible alterations of the gravitational field in a superconductor 
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FIG. 3: G. A. Ummarino, Possible alterations of the gravitational field in a superconductor 



